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ABSTRACT 

Experiments,  under  different  conditions,  were  conducted 
to  evaluate  the  effect  of  condensate  inundation  on  condensation 
heat  transfer  in  tube  bundles  of  marine  condensers.   Five  15.9 
mm.  (5/8  in)  nominal  outside  diameter,  smooth  stainless  steel 
tubes  were  used  in  a  vertical  row  to  simulate  an  actual 
condenser.   Tubes  were  located  in  an  equilateral  triangular 
array  with  a  spacing-to-diameter  ratio  of  1.5. 

Heat  transfer  performance  was  determined  for  each  tube  in 
a  bundle.   Data  was  taken  for  condensing  steam  on  the  outside 
of  each  tube  at  about  21  KPa  (3  psia)  and  at  about  101  PKa 
(14.7  psia).   Each  tube  was  cooled  by  water  on  the  inside  at 
velocities  of  0.78  to  6.22  m/sec  (2.56  ft/sec  to  20.42  ft/sec). 
The  overall  heat  transfer  coefficient  was  determined  directly 
from  experimental  data.   The  inside  and  outside  heat  transfer 
coefficients  were  determined  using  the  Wilson  plot  technique. 

Observation  of  condensate  flow  showed  lateral  droplet 
motion  along  the  first  three  tubes  in  portions  of  the  con- 
denser under  all  conditions  tested.   Side  drainage  occurred 
only  over  the  third  and  fourth  tubes  at  a  condensation 
pressure  of  about  21  KPa.   The  dominate  mode  of  the  flow  at 
101  KPa  condensation  pressure  was  gravitational  flow.   Outside 
heat  transfer  coefficients  were  higher  than  expected  under  all 


conditions  when  compared  to  the  Nusselt  theory.   The  reason 

for  this  is  possibly  due  to  secondary  vapor  flow.  Recommendations 

to  improve  validations  are  provided. 
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I.  INTRODUCTION 

A.   BACKGROUND  INFORMATION 

Recent  improvements  in  turbine  machinery  and  boiler  design 
have  brought  about  an  increase  in  horsepower  to  weight  ratio  of 
marine  propulsion  systems.   However  there  has  been  no  similar 
improvement,  in  practice,  in  the  steam  plant  condenser  size. 

Practical  marine  steam  condenser  design  is  based   almost 
exclusively  upon  two  documents.   The  Heat  Exchange  Institute 
(HEI)  standards  for  steam  surface  condensers  [1]  and  the 
standards  of  Tubular  Exchange  Manufacturers  Association 
(TEMA)  [2] .   These  standards  have  proven  to  be  reliable.   But 
it  is  evident  that  the  resulting  condensers  are  significantly 
overdesigned.   Briefly,  Search  [3]  has  shown  that  heat 
transfer  enhancement  methods  could  decrease  condenser   space 
to  weight  ratio,  thereby  establishing  new  design  criteria 
for  marine  condensers. 

It  has  been  the  objective  of  past  research  by  Eshleman  [4] 
to  investigate  the  effect  of  condensate  inundation  on  heat 
transfer  in  a  horizontal  tube  bundle.   But,  his  research 
concentrated  on  designing,  constructing  and  validating  the 
test  facility. 

Since  the  publication  of  Nusselt's  well-known  theoretical 
paper  [5]  on  film  condensation  many  theoretical  and 
experimental  studies  have  occurred. 
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The  evaluation  of  the  heat  transfer  coefficient  from 
the  vapor  condensing  on  the  outer  surface  of  horizontal 
tubes  is  based  mainly  on  Nusselt1 s  theoretical  formulas. 

Nusselt  performed  the  derivation  and  found  that  the 
average  heat  transfer  coefficient  n"  for  n  tubes  located 
below  one  another  is, 


h  =  0.725 
n 


pf(pf"pvJ  hfgkf  g 


ny^D   (T   -T  ) 
f  o    sv  w 


where   k_ 

r 


w 


pf 
pv 


Thermal  conductivity  of  film  (W/m-  C 
Average  tube  wall  temperature  (  C) 

Film  density  (kg/m  ) 

3 
Vapor  density  (kg/m  ) 

Dynamic  viscosity  of  film  (kg/m-sec) 


(  l  ) 


g  :  Acceleration  of  gravity  (m/sec  ) 


For  the  uppermost  tube,  equation  (1)  becomes 


h        =   0.725 
Nu 


Pf(Pf"Pv)    hfgkf   * 


U  JD    (T      -    T       ) 
f  o      sv       w 


n  h 


(    2    ) 


Hence,  the  relationship  between  the  average  heat  transfer 
coefficient  for  a  horizontal  tube  bundle  (h  )  consisting  of 
n  tubes  and  the  mean  heat  transfer  coefficient  (h^)  for  the 
uppermost  tube  may  be  found  as, 


5  /  h*T 

Nu 


-h 


(  3  ) 


The  above  equations  were  derived  by  Nusselt  assuming: 

1.   Condensation  of  a  saturated  vapor  at  negligible 
velocity, 
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2 .   Laminar  flow  of  the  condensate  film  in  a  continuous 
sheet  from  one  tube  to  "the  next  at  a  constant 
temperature  difference  for  all  tubes  in  the  bank. 

In  actual  condensers,  vapor  moves  at  a  fairly  high  speed 
over  a  considerable  part  of  its  path.   Under  changing  turbine 
speeds,  steam  velocity  is  not  negligible.   But  attempts  to 
evaluate  analytically  the  effect  of  vapor  velocity  have  not 
been  successful.   The  vapor  velocity  causes  friction  between 
the  vapor  and  the  condensate  film.   With  downward  flow  of 
vapor  and  condensate  (as  in  our  experiment)  the  frictional 
forces  are  added  to  the  force  of  gravity.   Consequently  as 
the  film  velocity  increases,  the  thickness  of  the  film 
decreases  and  the  coefficient  of  heat  transfer  from  vapor 
to  wall  increases . 

Additonally,  the  condensate  does  not  flow  down  from  each 
tube  in  a  continuous  film.   Rather  it  forms  drops  or  streams. 
Condensate  dripping  on  a  tube  from  above  splits  around  the 
tube  but  does  not  flow  axially.   The  thickness  of  the  film 
caused  by  the  condensate  coming  down  from  above  is  thus 
confined  to  the  place  where  the  drops  and  streams  descend, 
which  causes  local  disturbances  in  film  flow.   Briefly,  the 
true  nature  of  condensate  flow  on  and  between  tubes  differs 
from  the  flows  assumed  by  Nusselt. 

Recently,  Eissenberg  [6]  noted  that  the  condensate  droplets 
formed  on  the  tubes  strike  anywhere  on  the  lower  half  of  the 
tubes  below.   His  experimental  results  gave  heat  transfer 

coefficients  well  above  those  predicted  by  Nusselt. 
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Eissenberg  proposed  a  side  drainage  model,  and 
formulated  that, 

-h 
h  /  h   =  0.6  +  0.42  n  (  4  ) 

n    Nu 

As  a  matter  of  fact,  it  can  be  said  that  in  most  cases 
experimental  data  has  been  fit  to  various  modified  forms  of 
equation  (5)  below: 

h  /  h =  n"S  (  5  ) 

n    Nu 

where  0  .07£s<_0  .223  . 

Generally,  the  effects  of  condensate  inundation  and  vapor 
velocity  are  described  separately  from  each  other.   Actually, 
they  occur  simultaneously  and  their  combined  effect  is 
complex.   Fujii  [7]  correlated  the  effect  of  inundation  and 
vapor  shear  by  using  experimental  data  of  Nobbs  and  Mayhew 
[8,  9].   The  data  for  in-line  tube  banks  resulted  in 
equations  (6)  and  (7)., 

Nu°  =  10.74  Re?'312  (  6  ) 

m  L 

Equation  (6)  above  is  an  experimental  equation  for  a  tube 
without  inundation  where 

o 

Nu   :  Nusselt  number,  (h0dQ/kL)  for  pure  steam,  also 
without  inundation. 

Re-  :  Two  phase  Reynolds  number,  (U  d  /v.) 

kL  :  Thermal  conductivity  of  liquid 

d   :  Tube  outside  diameter 
o 

v.  :  Kinematic  viscosity  of  liquid 

and    U  :  Vapor  velocity 
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Nu  /Nu    is  comDUted  for  the  data  with  inundated  tubes 
mm 


and  equation  (7)  is  derived. 


Nu  /Nu°  =  (  Re/2xl06  ) 0 • 071 (w/w° ) ° ' 6 
mm        L 


(  7  ) 


Nu   :  Mean  Mussel t  number  for  a  tube 
m 

w   :  Rate  of  inundation  falling  onto  a  tube 
w   :  Rate  of  condensation  of  a  tube  correspond- 
ing to  Nu,  . 

3      m 

For  the  staggered  tube  bank,  a  correlation  was  not 
obtained  by  Fujii.Fujii  determined  that  for  the  staggered 
tube  bank,  the  inundation  effect  was  smaller  than  that  for 
the  in-line  bank. This  would  follow  from  examination  of 
Eissenberg ' s  side  drainage  model. 

Chisholm  [10]  combined  the  developments  of  Berman  and 
Tumanov  [11]  ,  and  Fuks  [12]  into  one  equation  to  evaluate  the 
effects  of  condensate  inundation,  vapor  velocity  and  non- 
condensable  gas  in  tube  bundles.   Chisholm  has  given  the 
following  formula  for  heat  flux, 


q    =     C.       (T        -     T     ) 
^  h  s  w 


3/4 


(     8     ) 


where 


C,     =    0.72  5 
h 


K       L    g 
c  3 


D    y     V 
o    c    c 


1/4 


n-1 


c  ,r 


1    + 


w 


-0.07 


>     (1  +  0.0095     ReU-8/VNu)     (9) 
v  m 


c  ,  n-1 


,/nu 
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Equation  (9)  is  for  the  downward  flow  of  the  vapor. 


w 


T 
s 

T 
w 

K 
c 

L 
yc 

V 

c 
c,r 

Ren, 
m 


Temperature  of  condensate  surface  (°C) 

Temperature  of  outer  tube  wall  (°C) 

Thermal  conductivity  of  condensate  (W/m-°C) 

Specific  latent  heat  (J/kg) 

Absolute  viscosity  of  condensate  (kg/m-sec) 

Specific  volume  of  condensate  (m  /kg) 

Rate  of  condensation  of  rth  tube  row  (kg/sec) 

Reynolds  number  of  vapor-gas  mixture. 


In  the  above  equation, 


Nu  = 


a    D 
stat  o 


K 


(  10  ) 


where  a  .  .  :  Heat  transfer  coefficient  across  the  condensate 
stat 

film,  static  vapor  condition  (W/m  -  C) 

Experimental  data  resulting  from  varying  inundation 
rates  and  vapor  velocities  is  sufficiently  scattered  to 
suggest  that  no  existing  correlation  fits  all  the  available 
data  (some  of  the  available  data  is  related  to  enhanced 
tubes) .   This  is  because  there  are  so  many  variables  that 
affect  marine  condenser  performance,  some  of  which  have  to 
be  discovered  (e.g.,  non-condensable  gas  effects,  pitch-to- 
diameter  ratio  of  condenser  tubes,  direction  of  vapor  flow, 
and  the  effects  due  to  ship  motion  in  three  dimension) . 
Research  is  presently  being  conducted  in  the  United  Kingdom 
aimed  at  producing  a  model  of  condensation  suitable  for 
sophisticated  condenser  performance  calculations.   However, 
the  results  have  yet  to  be  published. 
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B.   OBJECTIVES  OF  THIS  WORK 

In  order  to  evaluate  the  effect  of  inundation  on  a  bank 
of  enhanced  tubes,  it  was  necessary  to  establish  data  on  a 
bank  of  smooth  tubes  as  a  standard  of  comparison.   The 
objectives  of  this  work  were  therefore  :  (1)  to  establish 
experimental  data  on  a  bank  of  smooth  tubes  as  a  standard, 
and  (2)  to  compare  the  established  data  to  theoretical 
predictions . 
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II.   EXPERIMENTAL  APPARATUS 

A.  INTRODUCTION 

The  existing  test  facility  was  designed  by  Beck  [11]  and 
built  and  tested  by  Pence  [12] .   Major  modifications  were  made 
to  the  original  apparatus  by  Eshleman  [4] .   During  this  work, 
some  minor  modifications  were  made  and  these  changes  are  given 
in  Tables  II  and  III  for  each  run. 

B.  STEAM  SYSTEM 

The  steam  system  is  shown  in  Figure  2.   The  supply  of  steam 
is  locally  generated  and  supplied  to  the  building  which  houses 
the  experimental  apparatus.   The  steam  is  provided  by  means  of 
a  19.05  mm.  diameter  line  and  a  steam  inlet  valve  (MS-2).   A 
compound  pressure  gage  is  located  just  prior  to  the  steam 
separator  which  monitors  the  supply  pressure  as  it  is  adjusted 
by  (MS-2) .   Following  the  steam  separator,  a  line  strainer 
provides  additional  protection  from  contamination.   After  the 
strainer,  the  steam  proceeds  through  a  31.75  mm.  diameter  line 
which  provides  for  two  possible  steam  paths.   The  primary  path 
for  system  operation  is  via  the  throttling  valve  (MS-3), 
through  a  desuperheater  and  into  the  test  condenser.   Inside 
the  condenser,  the  steam  is  condensed  on  the  test  tubes.   The 
steam  which  is  not  condensed  proceeds  via  the  vapor  outlet  on 
the  test  condenser  to  the  secondary  condenser.   The  secondary 
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steam  flow  path  is  used  to  accomplish  system  stabilization 
during  startup  and  to  control  the  mass  flow  rate  of  steam  to 
the  condenser  during  operation.   Steam  proceeds  via  (MS-4) 
directly  to  the  secondary  condenser.   All  steam  lines  except 
the  primary  path  downstream  of  (MS-3)  were  insulated  with  25.4  mm, 
thick  fiberglass  insulation. 

C.   TEST  CONDENSER 

The  test  condenser  is  shown  from  various  views  in  Figures  3 
and  4.   Steam  enters  via  the  top  and  proceeds  over  the  baffle 
separators  and  through  a  flow  straightener ,  which  is  covered 
with  three  layers  of  150  mesh  screen,  to  the  tube  bundle.   The 
condensate  collects  at  the  bottom  of  the  condenser  and  flows 
out  one  of  the  two  12.7  mm.  diameter  openings  at  either  end  of 
the  condenser  to  the  hotwell  where  it  can  be  collected  and 
measured. 

Three  separate  viewing  windows  each  203  mm.  by  140  mm.  by 
17.7  mm.  and  made  of  pyrex  plate  glass  had  been  installed  to 
provide  maximum  viewing  of  the  active  tubes. 

The  tube  sheet  arrangement  is  shown  in  Figure  4.   The  tubes 
were  arranged  in  a  typical  naval  condenser  spacing-to-diameter 
ratio  of  1.5.   They  were  15.9  mm.  OD,  1.14  mm.  thick,  304 
stainless  steel  tubes  that  had  cooling  water  passing  through 
them.   Although  typical  naval  condenser  tubes  are  made  of  90-10 
copper-nickel,  the  choice  of  304  stainless  steel  was  based  on 
"on  hand"  stock  and  the  fact  that  the  principles  of  inundation 
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do  not  depend  on  the  tube  material  although  perhaps  the  heat 
flux  may  change  due  to  different  wall  resistance.   The  remaining 
half  tubes  were  made  of  15.9  mm.  OD  aluminum  bar  stock  and  were 
fastened  by  screws  to  the  outside  wall  of  the  steam  flow  path. 
This  arrangement  was  selected  to  best  simulate  the  steam  flow 
conditions  in  a  section  of  an  actual  condenser.   The  five  test 
tubes  are  singularly  removable.   The  top  tube  can  be  replaced 
by  a  304  stainless  steel  porous  tube  which  could  simulate  various 
condensate  inundation  rates. 

The  test  condenser  was  insulated  with  a  25.4  mm.  thick  sheet 
of  Armorflex  insulation. 

D.   CONDENSATE  SYSTEM 

The  condensate  system  is  shown  in  Figure  5.   The  condenser 
hotwell  collects  the  condensate  from  the  test  tubes,  while  the 
secondary  condenser  hotwell  collects  the  condensate  from  the 
secondary  condenser.   Valve  (C-l)  allows  the  isolation  of  the 
test  condenser  hotwell  for  condensate  measurement.   Since  house 
steam  was  used  as  the  steam  supply  system,  the  condensate 
collected  in  the  hotwells  is  pumped  back  to  the  house  system 
by  the  condensate  pump  via  valve  (C-3) .   The  condensate  lines 
were  insulated  with  19.1  mm.  thick  Johns-Manville  Aerotube 
insulation. 
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E.  COOLING  WATER  SYSTEM 

The  cooling  water  system  is  shown  in  Figure  6.   The  water 
used  was  normal  house  water  which  had  been  passed  through  a 
water  softener  on  the  way  to  the  supply  tank.   The  water  is 
pumped  from  the  supply  tank  by  a  5  HP  electric  driven  pump. 
It  is  routed  to  the  flowmeter  header  via  51  mm.  OD  plastic 
pipe.   The  flow  of  cooling  water  for  each  test  tube  is  then 
individually  controlled  and  measured  by  it's  own  rotometer. 
Each  rotometer  allowed  a  maximum  flow  rate  of  70.4  LPM.   The 
heated  cooling  water,  after  passing  through  the  test  section, 
was  piped  back  to  a  supply  tank.   A  separate  system  pumped  this 
heated  water  through  a  filter  and  cooling  tower  returning  the 
cooled  water  to  the  supply  tank  in  an  effort  to  maintain  a 
constant  cooling  water  inlet  temperature. 

After  leaving  the  rotometers,  the  system  piping  was  reduced 
to  15.9  mm.,  ensuring  a  distance  of  1  m.  ahead  of  the  test 
section  to  ensure  a  hydrodynamically  fully  developed  velocity 
profile  while  passing  through  the  test  section. 

F.  SECONDARY  SYSTEMS 
1.   Vacuum  System 

The  vacuum  in  the  test  condenser  and  secondary  condenser 
was  maintained  by  a  mechanical  vacuum  pump  and  a  vacuum  regulator 
which  induces  air  into  the  system.   The  vacuum  pump  takes  a 
suction  from  the  secondary  condenser  hotwell  which  is  connected 
to  the  test  condenser  hotwell  via  discharge  piping.   A  cold 
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trap  at  the  inlet  of  the  vacuum  pump  forces  incoming  vapor 
to  pass  over  a  system  of  refrigerated  copper  coils.   This 
removes  steam  and  entrained  water  from  the  vacuum  line  prevent- 
ing moisture  contamination  of  the  vacuum  pump  oil.   The  vacuum 
pump  outlet  is  vented  through  a  roof  exhaust  fan  to  avoid  a 
health  hazard  from  breathing  any  oil  vapor  exhausted  by  the 
pump.   A  schematic  diagram  of  this  system  can  be  found  in 
Figure  7. 

2 .   Desuperheater  System 

The  desuperheater  removes  sensible  heat  from  the  super- 
heated steam  by  injecting  water  at  about  25 °C  via  the  existing 
feedwater  system  through  valve  (DS-1)  and  a  rotometer.   The 
desuperheater  is  a  267  mm.  diameter  stainless  steel  can,  457  mm. 
high,  having  four  nozzles  inserted  equidistant  around  the 
circumference  of  the  inner  top  of  the  can.   The  nozzles  are  a 
fan  type  and  are  positioned  such  that  the  spray  is  downward  at 
a  45°  angle  to  allow  for  better  mixing.   A  collection  tank  is 
located  on  the  bottom  of  the  desuperheater  to  allow  for  drainage 
of  condensate.   This  system  can  be  isolated  by  valve  (DS-2) 

G.   INSTRUMENTATION 
1 .   Flow  rate 

a.   Steam  velocity  was  determined  by  calculation 


m  v 
c 


u  = 


s  <  X1  > 

A 
c 
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where 

m  '  =  Mass  flow  rate  of  condensate  (kg/sec) 

=  Q/hfg 

2 
A  =  Cross  sectional  flow  area  (m  ) 
c 

3 
and   v  =  Specific  volume  of  vapor  (m  /kg) 

b.   Cooling  water  flow  rate  was  measured  individually 

for  the  five  separate  tubes.   Each  flow  rate  was  determined 

by  a  rotometer  with  a  capacity  of  70.4  LPM  (18.6  GPM) .   The 

calibration  procedure  used  was  identical  to  that  listed  in 

Appendix  A  of  Ref .  [14] . 

2 .  Pressure 

Two  different  pressure  sensing  devices  were  used  during 
experimentation.   They  were  a  Bourdon  tube  pressure  gauge  which 
measured  steam  pressure  and  an  absolute  pressure  transducer 
coupled  with  a  760  mm.  mercury  manometer  which  was  used  to 
measure  test  condenser  pressure. 

3 .  Temperature 

Stainless  steel  sheathed  ,  copper-constantan  thermo- 
couples were  used  as  the  primary  temperature  monitoring  devices, 
Figure  3  shows  the  location  of  five  vapor  thermocouples.   The 
remaining  30  themocouples  of  this  type  were  located  as  shown  in 
Figure  6,  six  on  each  tube,  two  measuring  cooling  water  inlet 
temperature  and  four  measuring  water  outlet  temperature. 
Calibration  procedures  of  the  thermocouples  were  identical  to 
those  listed  in  Appendix  A  of  Ref.  [16]. 
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4.   Data  Collection  and  Display 

An  Autodata  collection  system  was  utilized  to  record 
and  display  the  temperatures  in  degrees  Celsius  obtained  from 
the  primary  stainless  steel  thermocouples.   Table  I  lists  the 
channel  numbers  and  location  of  these  devices. 
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III.   EXPERIMENTAL  PROCEDURES 

A.   OPERATING  PROCEDURES 

1.  Preparation  of  condenser  tubes 

Prior  to  any  run,  each  tube  was  properly  prepared  to 
ensure  filmwise  condensation.   The  cleaning  procedure  for 
stainless  steel  tubes  is  listed  in  Appendix  A. 

2 .  System  operation  and  steady  state  conditions 

The  basic  operating  instructions  developed  by  Pence 
[14]  and  modified  by  Reilly  [15]  were  used.   The  only  differ- 
ence in  the  procedure  as  listed  in  Appendix  B  of  Ref .  16  was 
that  instead  of  one  cooling  water  flowmeter  to  adjust,  the 
experimenter  had  five  to  set  as  desired. 

In  general  it  took  about  three  hours  from  initial 
light-off  until  steady-state  conditions  were  established. 
The  parameters  used  in  determining  steady-state  conditions 
were  cooling  water  inlet  temperature  and  steam  inlet 
temperature.   If  the  cooling  water  inlet  did  not  vary  more 
than  ±0.6  °C/hr  and  the  steam  temperature  did  not  vary  more 
than  ±0.3  °C/min,  steady-state  was  considered  achieved. 

The  time  for  the  system  to  stabilize  between  changes 
in  cooling  water  flow  rate  during  the  Wilson  plot  technique 
was  approximately  twenty  minutes.   This  time  increment  is 
suspect  as  other  investigators  waited  about  one  hour  for 
stabilization  between  changes,  especially  for  atmospheric  runs 
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It  must  be  pointed  out  however  that  the  amount  of  time  required 
to  collect  data  over  five  tubes  in  a  system  that  can't  be  shut 
down  and  repeated  the  next  day  prohibits  the  greater  time 
increment  between  data  points  for  the  Wilson  plot. 

The  general  set  up  for  the  data  taken  in  this 
research  was  a  steam  velocity  of  between  approximately  1.2 
m/sec  and  1.4  m/sec,  steam  temperature  of  62  °C  or  75  °C  for 
2lKPa  condensation  conditions  and  100  °C  for  atmospheric 
(101  KPa)  runs. 

B.   DATA  REDUCTION  PROCEDURE 

The  raw  data  collected  for  each  tube  for  run  10  can  be 
found  in  the  Tables  beginning  on  page  38. 

Appendix  B,  the  sample  calculations,  is  a  complete  listing 
of  the  equations  used  to  evaluate  the  data.   Appendix  C  is 
a  derivation  of  the  probable  error  in  the  data  reduction 
equations,  followed  by  a  sample  error  analysis  for  tube 
number  1  at  40  percent  flow,  run  10. 

The  following  standard  heat  transfer  equations  were  used 
to  reduce  the  raw  data  into  a  form  that  can  be  used  for 
evaluation. 

1.   Overall  heat  transfer  coefficient  (Un) 

The  method  employed  to  arrive  at  the  overall  heat 
transfer  coefficient  is  straightforward  and  similar  to  that 
employed  by  many  researchers  in  the  past.   The  heat  transfer 
rate  to  the  cooling  water  is  given  by 
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Q  =  m  C   (  T    -  T  .  )  (12) 

p    CO    ci  v  xz  ; 

The  heat  transfer  rate  from  the  steam  is  given  by, 

Q  =  m      C   (T  -T    )    +    h   +  C  ,    , (T    -  T    )      (13) 

con  |_pv  v  sat     fg  p(con)v  sat   con  J    K      ' 

The  heat  transfer  rate  can  also  be  found  from  the  overall 
heat  transfer  coefficient  by 


Q  =  Un  An  (LMTD)  (  14  ) 


where 


(t  -  t  . )  -  (t  -  t  ) 

V     Cl        V     CO 

LMTD  = (  15  ) 

Ln   ( T  -  T  .  )  /  ( T  -  T   ) 
v    ci     v    co 

After  combining  equations  (12),  (14),  and  (15)  it  can  be 
found  that 

T   -  T 
m  C  v    ci 

U  =  — —£ In  (  16  ) 

n      A  T   -  T  (  ' 

n  v    co 

2.   Corrected  overall  heat  transfer  coefficient  (Ur) 


u  = 


l 


c     1 


(  17  ) 


-  R 
U       w 

n 


where  R   is  the  wall  resistance  corresponding  to  different 
w 

tube  materials  and  may  be  given  by 

A  ln(r  /r. ) 
n    o   i 

Rw  =  (  13  ) 

2  IT  k  L 

w  ts 

where  L.   is  the  length  of  the  tube. 

Equation  (17)  allows  for  the  comparison  of  tubes  of  different 
materials  for  the  same  steam  and  cooling  water  condition 
within  the  test  condenser. 
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3.   Inside  heat  transfer  coefficient  (hi) 

h.  D. 
Nu  =  -U-  =0.036  Re0'8  Pr1/3  (L/D)"0'054  (  19  ) 


Equation  (19)  was  selected  because  both  the  Dittus-Boelter 
and  Sieder-Tate  relationships  which  are  commonly  used 
assume  a  fully  developed  velocity,  as  well  as  a  fully  devel- 
oped thermal  profile.   In  this  research,  it  was  suspected 
that,  although  the  velocity  profile  was  believed  to  be  fully 
developed,  the  thermal  profile  was  not  fully  developed. 
When  an  L/D  ratio  of  57.6  is  used  in  equaiton  (19)  a  constant 
of  0.029  results.   This  was  validated  by  computing  the  aver- 
age of  all  the  tube  constants  obtained  as  a  result  of  the 
Wilson  plot  technique.   Wilson  plots  for  each  tube  for  run  10 
can  be  found  in  Figures  8  through  12. 

4.   Outside  heat  transfer  coefficient  (hQ) 

The  outside  heat  transfer  coefficient  is  the  parameter 
that  is  used  to  compare  results  of  each  tube  in  the  bundle  and 

is  given  by 

hQ  =  (  20  ) 

(1/U  )  -  R   -  (D  /D.h.) 
n     w     oil 

Two  very  important  assumptions  were  made  in  using  this 
equation. 

a.   The  resistance  due  to  fouling  was  equal  to  zero. 
This  assumption  is  supported  by  the  fact  that  the  tubes  were 
new,  chemically  cleaned  and  smooth.   Also,  treated  soft  water 
was  used  as  the  cooling  medium. 
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b.   The  resistance  due  to  non-condensable  gases  was 
equal  to  zero.   This  assumption  was  supported  by  the  fact 
that  the  system  was  tested  for  air-tightness  and  found  to  be 
secure.   In  addition,  it  was  believed  that  the  velocity  of 
steam  passing  through  the  test  section  was  sufficiently  large 
to  keep  the  system  purged  of  any  non-condensables  that  might 
collect  in  the  test  section. 

C.   DATA  REDUCTION  COMPUTER  PROGRAM 

Reilly  [15]  developed  the  existing  program  in  Fortran 
Language.   His  program  had  been  translated  into  Basic  Langu- 
age for  use  with  the  HP  984  5  computer  by  Eshleman  [4]  during 
his  work.   Ultimately,  this  program  with  minor  modifications, 
can  be  used  in  an  integrated  system  between  the  Autodata  Nine 
data  collector  and  the  HP  984  5  computer.   This  will  allow 
automatic  data  input  with  real  time  data  output  for  the 
experimenter.   During  this  work,  the  existing  computer 
program  of  Eshleman  [4]  for  reduction  of  data  was  used  with 
the  HP  984  5  computer. 
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IV.  RESULTS  AND  DISCUSSION 

The  experiments  were  done  by  using  two  different  conden- 
sation pressures  to  establish  experimental  data  on  a  bank  of 
smooth  tubes .Condensation  pressure  for  runs  1  through  5  and 
9  through  12  was  maintained  at  about  3  psia.   For  runs  6,  7 , 
and  8  the  pressure  was  maintained  at  atmospheric  conditions. 
Runs  8,  9,  and  2  were  repetitions  of  runs  7,  5,  and  1 
respectively .Experimental  conditions  are  given  in  Tables  II 
and  III  for  each  run. 

In  equation  (16) ,  the  T   term  stands  for  vapor  temperature 
Large  differences  in  the  heat  transfer  coefficients  were 
obtained  depending  on  whether  the  actual  vapor  temperature  or 
the  saturation  temperature  was  used  in  equation  (16) .   It  is 
worthwhile  to  note  that  the  outside  heat  transfer  coefficients 
using  the  actual,  superheated  steam  temperatures  were  52  per- 
cent lower  than  those  calculated  using  the  saturation 
temperature. The  percent  change  of  U   values  was  7  9  percent. 
To  overcome  this  interesting  result,  it  was  decided  that  the 
T   term  in  equation  (16)  must  be  the  saturation  temperature 
instead  of  actual  vapor  temperature.  For  comparison  U   and 
hn  values  are  tabulated  in  Table  XIV  for  runs  10  (saturated) 
and  11  (superheat)  in  which  T  was  actual  vapor  temperature. 
Saturation  temperature  was  therefore  used  in  equation  (16) 
for  all  runs  except  those  reported  in  Table  XIV. 
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The  runs,  conducted  at  3  psia  condensation  pressure,  gave 

unexpected  results.   Outside  heat  transfer  coefficients  for 

the  five  test  tubes  showed  significant  fluctuations.   The 

values  of  h  decreased  for  the  first  three  tubes  and  then 
n 

increased.   Generally,  for  all  runs  conducted  at  3  psia,  the 
heat  transfer  coefficients  followed  this  same  pattern.   The 
cause  for  this  behavior  may  be  due  to  several  phenomena. 

During  observation  of  the  condensate  flow  pattern  on  the 
first  three  tubes,  there  was  evidence  of  lateral  droplet 
migration.   This  migration  was  due  presumably  to  a  deflection 
of  these  tubes  or  due  to  axial  flow  of  vapor.   This  causes 
a  non-uniform  heat  transfer  rate  across  the  length  of  the 
tube  because  of  a  decrease  in  film  thickness  at  some  locations 
and  then  an  increase  in  film  thickness  in  other  locations. 
This  thickening  of  the  film  on  the  tubes  would  result  in  a 
lower  average  heat  transfer  coefficient  than  expected. 
Secondly,  steam  appeared  to  concentrate  at  the  bottom  of  the 
test  condenser,  and  may  have  caused  cross  flow  around  the 
lower  tubes.   This  may  have  been  caused  by  a  system  resist- 
ance to  the  flow  of  excess  steam  which  is  not  condensed  in 
the  test  condenser.   Circulating  flow  of  the  excess  steam 
which  could  not  easily  leave  the  test  condenser  may  have 
caused  side  drainage  of  the  condensate.   It  was  observed 
that  the  condensate  flow  path  on  the  fourth  and  fifth  tubes 
was  toward  the  observation  window.   Because  of  these  two 
reasons,  the  outside  heat  transfer  coefficients  for  the  last 
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two  tubes  would  be  high.   As  an  example,  Table  XXII  lists 
values  for  run  10,  and  these  are  plotted  in  Figure  13  along 
with  two  theoretical  results. 

Figures  8  through  12  are  the  Wilson  plots  that  assist  in 
determining  the  constant  in  equation  (19) .   As  an  example, 
run  10  was  chosen  and  the  results  are  plotted  in  these 
Figures.   All  tubes  yield  good  linear  plots  with  slopes  which 
provide  constants  of  0.028,  0.032,  0.032,  0.028  and  0.030 
respectively.   The  expected  linear  plots  were  obtained  for 
all  the  other  runs  as  well  as  run  10.   The  data  reduction 
program  gives  the  option  of  using  the  constant  solved  for 
via  the  Wilson  plot  technique  or  inputing  one  of  the  user's 
own  choosing.   In  this  work,  for  all  runs,  0.029  was  used  as 
input  for  all  the  tubes  to  determine  the  inside  heat  transfer 
coefficient  which  in  turn  was  then  used  in  the  determination 
of  the  outside  heat  transfer  coefficient.   Outside  heat 
transfer  coefficients  for  all  runs  are  tabulated  in  Tables  XV 
through  XVIII. 

The  ratio  of  h  /h-,  as  listed  in  Tables  XIX  through  XXII 
for  all  runs,  was  determined  by  taking  the  average  outside 
heat  transfer  coefficient  h   for  n  tubes  and  averaging  them, 
then  dividing  by  the  outside  heat  transfer  coefficient  of 
the  first  tube  (h,  )  .   The  results  of  run  7,  as  listed  in 
Table  XXI,  are  plotted  on  Figure  14  along  with  the  theoretical 
equations  of  Nusselt  and  Eissenberg.   Based  on  the  observa- 
tions of  condensate  flow  at  atmospheric  pressure,  the  data 
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for  all  tubes  was  expected  to  fall  closer  to  the  Nusselt 
curve  due  to  the  presence  of  gravity  dominated  flow,  but 
certainly  not  below  it. 

The  experimental  study  on  the  effect  of  the  vapor  velocity 
upon  condensation  was  performed  both  at  3  psia  and  14.7  psia 
condensation  pressure.   Runs  9  and  10  at  3  psia,  and  runs  6 
and  7  at  atmospheric  pressure  were  conducted  at  different 
vapor  velocity  to  determine  the  effect  of  vapor  velocity. 
Vapor  velocities  for  runs  10  and  7  were  slightly  higher  than 
runs  9  and  6  respectively.   The  experimental  result  for  the 
3  psia  pressure  case  was  unexpected.   This  result  may  be  due 
to  the  effect  of  saturated  steam.   At  atmospheric  pressure, 
for  the  higher  vapor  velocity,  the  heat  transfer  coefficient 
was  higher  as  expected. 

From  examination  of  the  results  of  runs  7,  8  and  10,  11 
it  may  be  said  that  at  higher  condensation  pressure,  the  heat 
transfer  coefficient  is  higher. 

Each  repeated  run  was  within  the  uncertainty  range  of  the 
original.   The  conclusion  is  that  the  data  obtained  for  ten 
runs  is  sufficient  to  make  comparison  between  results. 

Side  pieces  and  baffles  were  installed  inside  the  test 
condenser  to  prevent  axial  flow  of  the  vapor.   However  droplet 
migration  still  resulted  in  both  cases.   From  the  comparison 
of  the  run  11  and  12,  it  may  be  concluded  that  side  pieces 
generally  did  not  have  any  effect  on  the  heat  transfer 
coefficient.   Heat  transfer  coefficients  of  runs  11  and  12 
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are  within  the  uncertainty  range  of  each  other.   The  differ- 
ence between  heat  transfer  coefficients  of  runs  4  and  10  is 
due  to  effect  of  baffles.   Outside  heat  transfer  coefficients 
are  higher  at  run  4  because  of  the  baffles  that  form  a  path 
for  the  condensate. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 
The  experimental  data  obtained  lead  to  the  following 
conclusions : 

1.  There  is  evidence  of  secondary  flow  of  steam  within 

the  test  condenser  which  is  suspected  to  have  influenced 
the  results. 

2.  Equation  (16)  should  not  be  used  with  superheated  steam 
temperatures . 

3.  The  measured  ratio  of  h n/h,  is  lower  than  expected  at 
atmospheric  runs. 

The  following  recommendations  are  provided: 

1.  Improve  the  steam  flow  path  to  ensure  a  uniform  downward 
profile  through  the  tube  bundle.   This  can  be  accomplished 
by  re-design  of  the  steam  inlet  section. 

2.  Prevent  auxiliary  system  resistance  against  the  flow 
of  uncondensed  steam  when  working  at  vacuum.   This  can 
be  done  by  either  conducting  experiments  at  atmospheric 
pressure  or  by  re-designing  the  auxilary  system  piping 
with  large  diameter. 

3.  Instrument  the  top  tube  with  thermocouples  to  measure 

T   in  order  to  check  the  outside  head  transfer  coefficient 
w 

against  the  Nusselt  theory. 

4.  Measure  the  vapor  pressure  between  tubes  to  get  more 
accurate  calculations. 
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TABLES 


TABLE    I 
LOCATION    OF    STAINLESS    STEEL    SHEATHED    COPPER-CONSTANTAN   THERMOCOUPLES 

CHANNEL  NUMBER  LOCATION  CHANNEL  NUMBER  LOCATION 


62 

64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 


hw 


shw 


v 


V 


T   . 

Cl 

§ 

1 

T   . 

Cl 

# 

1 

T   . 

Cl 

§ 

2 

T   . 

Cl 

# 

2 

T   . 

Cl 

§ 

3 

rn 

Cl 

§ 

3 

T 

CO 

§ 

1 

T 

CO 

# 

1 

T 

CO 

# 

1 

T 

CO 

# 

1 

T   . 

Cl 

# 

4 

T   . 

Cl 

§ 

4 

82 

T       #2 

CO 

83 

T       #2 

CO 

84 

T       #   2 

CO 

85 

T       #2 

CO 

86 

T       #   3 

CO 

87 

T       #3 

CO 

88 

T       #   3 

CO 

89 

T       #3 

CO 

90 

T   .    #  5 

Cl 

91 

T   .    #  5 

Cl 

92 

T       #4 

CO 

93 

T       #4 

CO 

94 

T       ^4 

CO 

95 

T       #4 

CO 

96 

T       §  5 

CO 

97 

T       #  5 

CO 

98 

T       US 

CO- 

99 

T       #  5 

CO 
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TABLE  II 
RUN  CONDITIONS  FOR  RUN  1  THROUGH  7 


Run  Mo . 

CONDITIONS 

1 

Condensation  pressure  3  psia ., Superheat 
steam, Low  vapor  velocity. 

2 

Same  as  Run  No:  1 

3 

Condensation  pressure  3  psia ., Superheat 
steam, Low  vapor  velocity , Steam  gage 
has  been  changed. 

4 

Condensation  pressure  3  psia ., Saturated 

steam, Low  vapor  velocity , Baffles   and 

Vacuum  control  valve  have  been  installed, 
Throttle  valve  has  been  changed. 

5 

Same  as  Run  No:  M- 

6 

Condensation  pressure  atmospheric, 
Saturated  steam, Low  vapor  velocity, 
Baffles  have  been  removed. 

7 

Condensation  pressure  atmospheric, 
Saturated  steam, High  vapor  velocity. 
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TABLE  III 
RUN  CONDITIONS  FOR  RUN  8  THROUGH  12 


Run  Mo . 

CONDITIONS 

8 

Same  as  Run  No .  7 

9 

Condensation  pressure  3  psia , Saturated 
steam, Low  vapor  velocity. 

10 

Condensation  pressure  3  psia. , Saturated 
steam, High  vapor  velocity. 

11 

Condensation  pressure  3  psia ., Superheat 
steam, High  vapor  velocity. 

12 

Condensation  pressure  3  psia ., Superheat 
steam, High  vapor  velocity , Side  pieces 
have  been  installed. 
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TABLE  IV 


RAW  DATA  FOR  TUBE  NO:  1,  RUN  10 


%  FLOW 

T   .(°C) 

CI 

24.950 

T     (°C) 

CO 

36.067 

T     (°C) 

V 

73.375 

GPM 

10 

1.92 

12.5 

25.250 

34.325 

73.375 

2.40 

15 

25.750 

34.767 

73.375 

2.38 

17.5 

26.550 

34.767 

73.375 

3.36 

20 

27.100 

34.467 

73.375 

3.84 

25 

26.450 

32.200 

73.375 

4.80 

40 

27.700 

32.700 

73.375 

7.68 

50 

28.700 

33.167 

73.375 

9.67 

70 

30.100 

33.533 

73.375 

13.48 

80 

31.050 

34.133 

73.375 

15.36 
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TABLE  V 


RESULTS  FOR  TUBE  NO:  1,  RUN  10 


U 


n 


U 


Q 


m 


2906.755 
2912.181 
3509.480 
3763.169 
3864.645 
3644.438 
5169.389 
5912.132 
6463.177 
6736.134 


3486.533 
3494.363 
4391.116 
4795.622 
4961.648 
4604.459 
7340.180 
8933.859 
10255.078 
10959.733 


5649.880 

5767.253 

6874.809 

7307.513 

7487.057 

7310.103 

10166.006 

11431.371 

12241.189 

12521.407 


.002405 
.002455 
.002926 
.00  3111 
.003187 
.00  3112 
.004327 
.004866 
.005211 
.005330 


h. 

i 

h 
o 

23811.652 

C 

Re 

C 
wv 

4598.834 

0.029 

14054.5 

.78 

5453.646 

12542.547 

0.029 

17314.4 

.97 

6343.280 

19904.566 

0.029 

20975.9 

1.17 

7207.769 

19113.540 

0.029 

24669.7 

1.36 

8031.485 

16296.772 

0.029 

28264.7 

1.55 

9446.257 

10204.903 

0.029 

34305.3 

1.94 

13893.319 

18116.447 

0.029 

55869.9 

3.11 

16841.966 

22180.994 

0.029 

71391.4 

3.91 

22183.212 

21386.404 

0.029 

101287.7 

5.46 

24834.745 

21782.715 

0.029 

117197.2 

6.22 
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TABLE  VI 


RAW  DATA  FOR  TUBE  NO:  2,  RUN  10 


%  FLOW 

T   .(°C) 
Cl 

25.000 

T     (°C) 

CO 

35.425 

T     (°C) 

V 

73.375 

GFM 

10 

1.92 

12.5 

25.300 

34.325 

73.375 

2.40 

15 

25.850 

34.150 

73.375 

2.88 

17.5 

26.600 

34.175 

73.375 

3.36 

20 

27.100 

33.825 

73.375 

3.84 

25 

26.250 

32.450 

73.375 

4.80 

40 

27.750 

32.025 

73.375 

7.68 

50 

28.800 

32.500 

73.375 

9.67 

70 

30.200 

32.850 

73.375 

13.48 

80 

31.100 

33.500 

73.375 

15.36 
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TABLE  VII 


RESULTS  FOR  TUBE  NO:  2,  RUN  10 


U 


n 


U 


Q 


m 


2705.496 
2897.657 
3209.670 
3446.0  36 
3500.473 
39  32.755 
4387.397 
4863.757 
4953.877 
5206.388 


3200.943 
3473.472 
3931.614 
4292.242 
4377.025 
5074.475 
5857.697 
5738.890 
6913.138 
7415.001 


5298.997 

5735.404 
6328.988 
6737.508 
6835.719 
7882.096 
8693.329 
9469.938 
9450.620 
9748.893 


.002256 
.002441 
.002694 
.002868 
.002910 
.003355 
.003700 
.004031 
.004023 
.004150 


Re 


C 


wv 


4583.724 

5455.170 

6325.033 

7186.119 

8002.946 

9448.909 

13844.978 

16789.189 

22112.360 

24756.032 


149  74.334 
12268.806 
13098.673 
13105.593 
11391.899 
12835.402 
11186.354 
12295.258 
10668.317 
11187.876 


0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 


13971.0 

.78 

17323.2 

.97 

20866.8 

1.17 

24535.6 

1.36 

28083.1 

1.55 

34  322.7 

1.94 

55518.4 

3.11 

70986.6 

3.91 

100702.3 

5.46 

116524.8 

6.22 
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TABLE  VIII 


RAW  DATA  FOR  TUBE  NO:  3,  RUN  10 


%  FLOW 

T   .(°C) 

CI 

25.025 

T     (°C) 

CO 

35.450 

T     (°C) 

V 

73.375 

Gm 

10 

1.92 

12.5 

25.300 

34.325 

73.375 

2.40 

15 

25.800 

34.125 

73.375. 

2.88 

17.5 

26.500 

34.100 

73.375 

3.35 

20 

27.100 

33.725 

73.375 

3.34 

25 

25.500 

32.400 

73.375 

4.80 

40 

27.750 

32.000 

73.375 

7.58 

50 

28.800 

32.500 

73.375 

9.67 

70 

30.200 

32.775 

73.375 

13.48 

80 

31.100 

33.425 

73.375 

15.36 

45 


u 


n 


TABLE  IX 


RESULTS  FOR  TUBE  NO:  3,  RUN  10 


U 


m 


2701.263 

3195.018 

5286.187 

.002250 

2897.657 

3473.472 

5735.404 

.002441 

3216.662 

3942.109 

6  348.174 

.002702 

340  8.819 

42  34.655 

6670.928 

.002840 

3444.271 

4289.504 

6734.245 

.002866 

3750.225 

4774.622 

7500.319 

.003193 

4360.472 

5809.799 

8642.547 

.00  3679 

4863.757 

6738.890 

9469.938 

.004031 

4809.36  3 

6634.918 

9183.326 

.003909 

5039.093 

7080.227 

9444.421 

.004020 

h. 

i 

h 
o 

4585. 

,638 

14822. 

,994 

5455. 

,170 

12268. 

,806 

6322. 

,386 

13228. 

,923 

7183. 

124 

1259  3. 

,467 

7998. 

,502 

10826. 

,610 

9459. 

,522 

11059, 

,603 

13843. 

,045 

11014. 

,343 

16789. 

,189 

12295, 

.258 

22103, 

,247 

10022, 

.034 

24745 

.908 

104444.89: 

Re 


C 


wv 


0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 


13981.6 

.78 

17323.2 

.97 

20851.0 

1.17 

24517.1 

1.36 

28054.8 

1.55 

34  392.6 

1.94 

55504.3 

3.11 

70986.6 

3.91 

100627.1 

5.46 

1164  38.4 

6.22 
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TABLE  X 


RAW  DATA  FOR  TUBE  NO:  4,  RUN  10 


%  FLOW 

T   .(°C) 
CI 

24.500 

T     (°C) 

CO 

35.200 

T     (°C) 

V 

73.375 

Gm 

10 

1.92 

12.5 

24.800 

34.067 

73.375 

2.40 

15 

25.300 

33.867 

73.375 

2.88 

17.5 

26.100 

33.833 

73.375 

3.36 

20 

26.600 

33.700 

73.375 

3.84 

25 

26.000 

33.100 

73.375 

4.80 

40 

27.200 

31.900 

73.375 

7.68 

50 

28.300 

32.367 

73.375 

9.67 

70 

29.700 

32.833 

73.375 

13.48 

80 

30.600 

33.533 

73.375 

15.36 
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u 


n 


TABLE  XI 


RESULTS  FOR  TUBE  NO:  4,  RUN  10 


U 


m 


2754.738 

3270.101 

5439.790 

.002  315 

2950.654 

3549.902 

5890.343 

.002507 

3282.571 

4041.559 

6533.981 

.002781 

3484.748 

4352.466 

6879.525 

.002928 

3670.281 

4645.789 

7218.049 

.003072 

4526.156 

6107.690 

9025.343 

.003842 

4787.990 

6594.309 

9559.236 

.004069 

5308.371 

7623.594 

10410.942 

.0044  32 

5322.347 

8730.418 

11174.606 

.004757 

6328.286 

9919.585 

11915.378 

.005072 

h. 

l 

h 
o 

4565.227 

16897.966 

5432.063 

13435.250 

6295.645 

14578.640 

7152.500 

13823.001 

7975.173 

13499.872 

9470.136     ■ 

22301.564 

13792.792 

14282.183 

16730.289 

15655.474 

22049.547 

15752.896 

24692.997 

18110.924 

Re 


C 


wv 


0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 


13869.1 

.78 

17190.5 

.97 

20691.5 

1.17 

24328.1 

1.36 

27906.8 

1.55 

34462.5 

1.94 

55140.0 

3.11 

70536.2 

3.91 

100184.6 

5.46 

1159  87.6 

6.22 
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TABLE  XII 


RAW  DATA  FOR  TUBE  NO:  5,  RUN  10 


%  FLOW 

T   .(°C) 
d 

24.900 

T     (°C) 

CO 

36.067 

T     (°C) 

V 

73.375 

GFM 

10 

1.92 

12. 

5 

25.150 

34.733 

73.375 

2.40 

15 

25.750 

34.500 

73.375 

2.88 

17. 

5 

26.450 

34.567 

73.375 

3.36 

20 

27.050 

34.267 

73.375 

3.84 

25 

26.400 

32.867 

73.375 

4.80 

40 

27.650 

32.433 

73.375 

7.68 

50 

28.700 

32.767 

73.375 

9.67 

70 

30.050 

33.167 

73.375 

13.48 

80 

31.000 

33.767 

73.375 

15.35 
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TABLE  XIII 


RESULTS  FOR  TUBE  NO:  5,  RUN  10 


U 


n 


U 


Q 


m 


2918.294 
3087.244 
3394.481 
3704.291 
3774.685 
4128.802 
4926.846 
5357.143 
5839.040 
6015.076 


3503.168 
3749.482 
4212.551 
4700.415 
4814.341 
5405.668 
6860.610 
7724.593 
8768.006 
9171.035 


5675.364 

6089.610 

6671.698 

7219.136 

7335.082 

8220.336 

9725.592 

10408.808 

11115.595 

112  39.183 


.002416 
.002592 
.002840 
.003073 
.003122 
.003499 
.004140 
.004431 
.004731 
.004784 


h 


Re 


wv 


4597.558 

5463.038 

6333.856 

7195.785 

8020.371 

9478.999 

13868.798 

16804.731 

22132.653 

24778.577 


24649.985 
16498.431 
16770.374 
17767.523 
14851.253 
15102.528 
15485.062 
16010.637 
15827.522 
15723.011 


0.029 

0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 
0.029 


14047.5 

.78 

17368.4 

.97 

20919.5 

1.17 

24595.4 

1.36 

2819  3.9 

1.55 

34520.9 

1.94 

55691.5 

3.11 

71105.7 

3.91 

100869.8 

5.46 

116717.2 

6.22 
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TABLE  XIV 


CALCULATED  U   AND  h  VALUES  FOR  RUN  10  AND  11  BY  USING  ACTUAL 
n      n 


VAPOR  TEMPERATURE 


Run  No:  10 
Saturated 

Run  No:  11 
Superheat 

U 
n 

h 
n 

U 
n 

h 
n 

4328 .7910 
3774.5253 
3715.2327 
4143.9355 
4163.9537 

15785.4551 
10963.0027 
10574.6464 
13969.3079 
14537  .9795 

3415. 5756 
2988 .8119 
2938  .9235 
3351.3064 
3290.4901 

7720.4862 
5985.2128 
5819.6943 
7418.9868 
7238  .0133 
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TABLE  XV 
OUTSIDE  HEAT  TRANSFER  COEFFICIENTS  FOR  RUN  1,  2  AND  3 


Run  No:  1   ' 

h  (W/m2-°C) 

n 

Run  No:  2 

h  (W/m2-°C) 
n 

Run  No :  3 

h  (W/m2-°C) 

n 

20410.5000 
9211.6182 
10419. 5087 
20683. 4820 
18038.4950 

22051.0505 
10891.4140 
11583 .6072 
19188 .1990 
16858.8026 

17982  .7509 
10026.5844 
10907.6969 
4420  .6795 
14247  .3599 

TABLE  XVI 


OUTSIDE  HEAD  TRANSFER  COEFFICIENTS  FOR  RUN  4,  5  AND  6 


Run  No:  4 

h  (W/m2-°C) 

n 

Run  No :  5 

h  (W/m2-°C)  " 

n 

Run  No :  6 

h  (W/m2-°C) 

n 

23968 .0185 
12797.6682 
13595.7268 
18680.0589 
19960.6512 

23357.4232 
11543 . 5844 
12902  .0927 
18936.9402 
19348 .2619 

21608  .6195 

11008 .3957 

10491.5847 

9614 .5084 

5027.7868 
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TABLE  XVII 
OUTSIDE  HEAT  TRANSFER  COEFFICIENTS  FOR  RUN  7 ,  8  AND  9 


Run  No:  7 

h  (W/m2-°C) 
n 

Run  Mo:  8 

h  (W/m2-°C) 
n 

Run  No:  9 

h  (W/m2-°C) 
n 

28588.3804 
13356.1272 
15115.2214  ' 
12932.0973 
12006.7075 

26265.9234 
14002 .0691 
12683 .1094 
11423.2917 
10005.0352 

21990.0288 
13282  .0589 
12600.4265 
16127.8443 
18357.4732 

TABLE  XVIII 
-OUTSIDE  HEAT  TRANSFER  COEFFICIENTS  FOR  RUN  10,  11  AND  12 


Run  Mo:  10 

h  (W/m2-°C) 
n 

Run  No:  11 

h  (W/m2-°C) 
n 

Run  No:  12 

h  (W/m2-°C) 
n 

18534.0540 
12301.2512 
11357 .6936 
15833.7770 
16868.6326 

16052.7204 
11114.5881 
10884 .2246 
14821.7962 
15112  .0203 

13965.0102 
10540.0521 
10139 .7836 
■13651.0698 
13800.1873 
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TABLE  XIX 


DATA  RESULTS  FOR  RUN  1,  2  AND  3 


Run  No:  1 

Run  Mo:  2 

Run  No :  3 

h  /h, 
n   1 

h  /h1 

n   1 

Vhi 

1.00000 

1.00000 

1.00000 

0.72566 

0.74696 

0  .77878 

0.65394 

0.67308 

0.72138 

0.74379 

0.72235 

0.60249 

0.77179 

0.73079 

0.64045 

TABLE  XX 

DATA  RESULTS  FOR  RUN  4,  5  AND  6 


Run  No :  4 

Run  Mo:  5 

Run  No :  6 

h  /h- 

n   1 

n   1 

h  /h, 

n   1 

1.00000 

1.00000 

1.00000 

0.76697 

0  .74711 

0.75472 

0.70040 

0.68220 

0.66499 

0.72014 

0.71433 

0.60998 

0.74267 

0.73714 

0.53452 
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TABLE  XXI 


DATA  RESULTS  FOR  RUN  7 ,  8  AND  9 


Run  No:  7 

h  /h- 

n   1 

Run  No:  8 

Vhi 

Run  No:  9 

h  /h, 

n   1 

1.00000 
0.82104 
0.72360 
0.65579 
0.60863 

1.00000 
0  .76654 
0.67199 
0.61272 
0.56636 

1.00000 
0.80200 
0.72567 
0.72761 
0.74  9.0  5 

TABLE  XXII 


DATA  RESULTS  FOR  RUN  10,  11  AND  12 


Run  No:  10 

Run  No:  11 

Run  No:  12 

h  /h, 
n   1 

h  /h, 

n   1 

h  /h, 

n   1 

1.00000 

1.00000 

1.00000 

0.83186 

0.84619 

0  .87737 

0.76783 

0 .79014 

0.82694 

0.78945 

0.82343 

0.36459 

0.81359 

0.84703 

0.83931 
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FIGURES 


Fig.  1.   Droplet  Path  Through  a  Tube  Bundle  with  Side 
Drainage 
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23.81mr 


7.94  mm 

V 


15.9  mm   00 
ACTIVE   TUBES    (STAINLESS   STEEL) 


DUMMY    TUBES   (ALUMINUM) 
15.9  mm    00 


SO    =1.5 


Fig.   4    Schematic  Side  View  of    Test   Tube 
Arrangement. 
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NUMBER        OF     TUBES CnJ 


Fig.  13.   Average  Outside  Heat  Transfer  Coefficient  Ratio 
Versus  Number  of  Tubes  for  Run  No: 10 
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m  co  ♦  in 

NUMBER      OF    TUBES    [n] 


'£ 


Fig.  14.   Average  Outside  Heat  Transfer  Coefficient  Ratio 
Versus  Numberof  Tubes  for  Run  No: 7 
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APPENDIX  A 
TUBE  CLEANING  PROCEDURE 

To  ensure  filmwise  condensation,  the  condenser  tubes  had 
to  be  prepared.   Surfaces  of  the  tubes  were  cleaned  to  insure 
proper  wetting  characteristics  and  to  insure  that  all  deposits 
were  removed.   Stainless  steel  tubes  were  prepared  in  accordance 
with  the  procedure  given  in  Newton  [17] .   The  steps  in  this 
cleaning  procedure  are  as  follows: 

1.  Prepare  an  Alconox  detergent  solution  and  heat  to  90  °C. 

2.  Apply  this  solution  to  the  surface  of  the  tubes. 

3.  Drain  and  rinse  the  tubes  with  distilled  water. 

4.  Spray  with  alcohol. 

5.  Rinse  with  distilled  water. 

6 .  Spray  with  acetone 

7.  Rinse  with  distilled  water. 
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APPENDIX  B 
SAMPLE  CALCULATIONS 
The  following  is  an  example  of  how  the  data  reduction 
program  calculates  the  results.   Tube  number  one  at  40  percent 
flow  rate  of  cooling  water  of  run  10  was  selected  for  this 
analysis.   This  same  tube  and  flow  rate  was  used  for  the  error 
analysis  in  Appendix  C. 

Input  parameters 

Tube  Outside  Diameter  (D  )  0.015875  m. 

o 

Tube  Inside  Diameter  (D.)  0.0141  m. 

Tube  Length  (L   )  0.9144  m. 

2 

Outside  Nominal  Surface  Area  (A  )  0.0456  m 

n 

Wall  Resistance  (R  )  5.72xl05  m2-K/W 

Cooling  Water  Inlet  Temperature  (T  . )  27.7  °C 

Cooling  Water  Outlet  Temperature  (T   )  32.7  °C 

Average  Cooling  Water  Temperature  (Tbc/Tbk)  30.2  °C,303.4  K 

Steam  saturation  temperature  (T  )  73.375  °C 

Gallons  Per  Minute  of  Cooling  Water  (GPM)  7.68  GPM 

Section  1.   Water  Properties 

y  (MHUW)  =  (4.134xl0~4)  Exp   |  [(0. 00829158)  (303 . 4) +  ( 2644. 2184) /( 303 .4)1 

-  10.59252566  > 

-4 
y        =  7.828  x  10     kg/m-sec 
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k(KW)  =  0.5565919  +  (0.002174417)  (30.2)  -  (0.70127  x  10~5)  (30. 2)  2 


-  (2.0914  x  10~10) (30. 2)3 


K  =  0.615858  W/m-°C 


P  (RHO)  =1004.44434-  (0 . 12673368)  (30 . 2 ) 

-  (  0.0023913147)  (30. 2) 2 


p    -     998.436   kg/m 


C        (CP)     =    4.2377955    -    (0.0018553514    )     (30.2) 

+     (1.3948314   x    10~5)(30.2)2 


C       =4.195   kJ/   kg-°C 


itl     (MFRCW)     =   LPM   x   BHC'   x   1'el   xl° 

where      LPM=GPM  x   3.78533 

jjj    =      (29.0713344)  (998.436)  (1.67xl0~5) 

m   _     0.484732  kg/sec. 

Prandtl   Number    (Pr) 

Pr    =     UC   /   k   =    (7.828   x   10~4x   4.1945x   103)    /    (0.615858) 
P 


5.3315 


Pr  = 


Section  2.  Data  Reduction 


2 
1.   Cooling  water  velocity  (C   )  =  4  m  /pir  d. 


C___  =   (4  x  0.484732)  /  j  (998.436)  TT  (0.0141)  2 


C         =     3.1092  m/sec 

CW 
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Mass  flow  rate  per  unit  area  (G) 
2 

(i    =4m/7TD.    =pC 
u  1  cw 

G    =   (998.436) (3.1092) 

2 
q    =  3104.3372  kg/m  -sec 

Reynolds  Number    (Re) 

Re  =D.G  /  ]i   =  (0.0141  x  3104.3373)/  (7.828  x  10  ') 
Re  _  55916.14 

Overall   Heat   transfer   Coefficient    (U    ) 

n 
m  c           ,     /  T     -  T    . 
U_    =  p         In  /      v ci_ 


-4, 


n 


T     -  T 

v  CO 


-) 


=      (0.484732) (4.1945    x    10    ) 
0.0456 

=     5169.408     W  /  m2    -   °C 


In 


(73.375    -    27.7 
73.375    -    32.7 


5.      Corrected   Overall   Heat   Transfer   Coefficient    (U    ) 

c 

•    U  1 1 

c  1 


u 


-      R 


W 


5169.408 


-    5.72    x   10 


-5 


U       =      7339.68  W   /  m2-    °C 

c 


6.      Wilson   Plot   Parameters    (X,Y) 
(a)    Ordinate 

Y   =  1 L_ 


u 

n 
y    =  19.345   x   10 

(b)    Abscissa 

X   = I 


-5 


X    = 


_     9.114   x   10 


-5 


5169.408 

2      o 

m  -      C   /  W 


_  0.8  „   1/3 

Re  Pr 


(55916. 14)°*8(5. 3315) 1/3 
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7.   Determination  of  constant 
D 

C  =    Mk~~ 

Where  M  =  Slope  returned  by  linear  regression  subroutine 

N  =  0.9265 

0.015875 


C  = 


(0.9265) (0.615858) 
C  =  0.028 


8.   Inside  Heat  Transfer  Coefficient  (h.) 

i 

Nu  =  -4^-  =  0.036  Re0'8  Pr1/3  (L/D  j"0'054 
k  o 

where 

L        0.9144 


0.015875 


=  57.6 


h.   =  — |—  (0.029)  Re0"8  Pr1/3 

l       D. 

l 

h.   =  °-615858    (Q_02    (55916.14)0.8         1/3 

i       0.0141 

h.   =  13898.64  w/m  -°C 
i 

9.   Outside  Heat  Transfer  Coefficient  (h  ) 

o 


h  = 
o 


h  = 
o 


D 

-  R   - 


U       w    D.h. 
n  ii 


0.015875 
-  5.72  x  10 


5169.408  (0.0141) (13898.64) 

h  =  18103.22  W/  m2-  °C 
o 
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APPENDIX  C 
ERROR  ANALYSIS 


The  basic  equations  used  in  this  section  are  reproduced 
from  Reilly  [15] .   The  general  form  of  the  Kline  and 
McClintock  [18]  "second  order"  equation  is  used  to  compute 
the  probable  error  in  the  results.   For  some  resultant,  R, 


which  is  a  function  of  primary  variables  x, ,  x~ 
the  probable  error  in  R,  6R  is  given  by: 


5R  = 


kn' 


->  U 


IT   ^ 
5x,   2 


6R  . 
5x2   2 


'  (*S 


,  5x   is  the  probable  error  in  each  of 


where  5x,,  5x2/  .... 
the  measured  variables. 

1.   Uncertainty  in  overall  heat  transfer  coefficient,  U 


n 


6U 


n 


U 


n 


SA  \  2    /  5C  \  2 


5m  \  2 


n 


m 


5T  (T  .-T   ) 
v   ci   co; 


-\2 


T  -T  . 
v  ci  , 


w (T  -T  • ) (T  -T   )  Ln  , 

v   ci    v   co      T..-T 


v   co 


6T  . 
ci 


,2 


+  T  -T 

_(T  -T  . )  Ln   v   CI  J 

v  cl    T-^r— 

V   CO 


5T 


co 


L  (T  -T   )  Ln  Tv"Tci  J 
v  co 


T  -T 

V   CO 


The  following  values  are  assigned  to  the  variables 
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5A   =  ±0.0001  m' 
n 


5C   =  ±0.0042  KJ/kg-°C 
5m  =  ±0.01  kg/sec 


5Tv  =  +0.5  °C 


6T  •  =  ±0.1  °C 
ci 


5T  n   =  ±0.1  °C 

CO 


For  tube  No:l  at  40  percent  of  10.  run 


5U 


n 


U 


n 


/p. oooi  N 

\^0.0456  J 


0.0042 
4.1945 


+ 


[ (0.5)  (-5) 1 

[(45.675)  (40.675)  Ln(1.1229)J 


0.01 

0.484732 


+ 


r        o.i        i2     r 

[  (40.675)  Ln(1.1229)        [   C45.67 


0.1 


1/2 


-»2 


5)  Ln  (1.1229),, 


5U 


n 


U 


=  0.037 


n 


2.   Uncertainty  in  inside  heat  transfer  coefficient,  h. 
The  probable  error  in  the  inside  heat  transfer 
coefficient  is  given  by: 
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6h 


5k 
k 


6Di  V  Vo.8  6Re 


Di 


Re 


0.333  6Pr 
Pr 


(¥) 


-■  U 


where 

6k  =  ±0.001  W/m-°C 


5D.  =  ±0.001  m 


<5Pr   =   ±0.10 
6C   =   +0.001 


6  Re 
Re 


6G 
G 


(0.01) A    + 


6y 


0.1 

7.828 


+ 


D. 

l 


2   "1 


0.001 
0.0141 


h 


5  Re 
Re 


=    0.073 


6h. 


0.001 
0.615858 


0.001 
0.0141 


+      0.333 


0.10   \  2      /o.ooi 
.3315y       +  \^0.029 


+    (0.8x0.073) 


=      0.098 


h.f    Tube   No:l_40%_10.Run=13899    ±   1362   W/m   -°C 
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3.   Uncertainty  in  outside  heat  transfer  coefficient,  h 
The  probable  error  in  the  outside  heat  transfer 
coefficient  is  given  by: 


6h 


=  < 


6Un 


n        D.h. 


+ 


5R 


w 


1  1 


-L  -  p  _   ° 

L   U     w   D.h.   J 
n         li 


D 


Sh. 

i 


Dihi 


_L  _  R 9_ 

L   u      w    D-h. 
n  ii 


where   6Un   =  0.037 


U 


n 


<SR  =  2.86x10  6  m2-°C/W 
w 


oh. 

l 


=  0.098 


U   "  Rw    D.h. 
n         ii 


=  5.524xl0"5  m2-°C/W 


6h. 


0.037 


(5169.408)  (5.524xl0-5) 


(8.10xl0~5)  (0.098) 


2.86x10 
5.524x10 


-6 

^5 


(5.524x10 


-5 
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5h 

T-2-  =  0.200 
h 
o 


h  ,  Tube  No:l-40%-10.Run  =  18103±3621  W/m2-°C 
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